Dewetting dynamics of stressed viscoelastic thin polymer films 
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Ultrathin polymer films that are produced e.g. by spin-coating are believed to be stressed since 
polymers are 'frozen in' into out-of-equilibrium configurations during this process. In the framework 
of a viscoelastic thin film model, we study the effects of lateral residual stresses on the dewetting 
dynamics of the film. The temporal evolution of the height profiles and the velocity profiles inside 
the film as well as the dissipation mechanisms are investigated in detail. Both the shape of the 
profiles and the importance of frictional dissipation vs. viscous dissipation inside the film are found 
to change in the course of dewetting. The interplay of the non-stationary profiles, the relaxing initial 
stress and changes in the dominance of the two dissipation mechanisms caused by nonlinear friction 
with the substrate is responsible for the rich behavior of the system. In particular, our analysis 
sheds new light on the occurrence of the unexpected maximum in the rim width obtained recently 
in experiments on PS-PDMS systems. 

PACS numbers: 68.60.-p,68.15.+e,83.10.-y 
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I. INTRODUCTION 

Thin polymer films arc not only of obvious techno- 
logical importance e.g. for coatings and lubrification 
purposes but also have attracted recently the attention 
of physicists because of their rich dynamical behavior 
[II, S 111- On a non-wettable substrate, thin polymer 
films are unstable and start to dewet. For usual, purely 
viscous liquid films this has been known and studied al- 
ready for some decades (see 0, Q for review articles on 
this subject), but for polymeric fluids the dynamics and 
phenomenology is by far richer. As thinner and thin- 
ner films became technologically feasible, films produced 
with thicknesses smaller than the equilibrium size of a 
single polymer have been studied [1,13, d, dl- These ex- 
periments revealed asymmetric rim shapes, several dis- 
tinct regimes for the temporal evolution of the dewetting 
velocity and a non-monotonous behavior of the width of 
the rim. 

The asymmetric rim shapes could be reproduced by 
models assuming either a shear thinning fluid [To'| or a 
viscoplastic solid 11]. Such nonlinearities in the mechan- 
ical properties might well be present, but seem to be less 
essential for the formation of the rim than viscoelasticity 
and the friction with the substrate [H, [H, [ij • This as- 
sertion is supported by the absence of rim formation in 
recent experiments on the dewetting of PS (polystyrene) 
films floating on a non-wettable liquid [l3| , where friction 
is avoided, while linear viscoelasticity could be clearly 
detected. In the experiments on supported films men- 
tioned above [1, [^, the silicium substrate is usually 
coated with PDMS (polydimethylsiloxane), leading to 
a polymer-polymer interface at the substrate inducing 
a strong slippage [l3|. Indeed, in a viscoelastic model 
assuming slippage and friction with the substrate, the 
asymmetric shapes could be reproduced [isj . 

Even more striking than the rim shapes has been the 
occurrence of a maximum in the width of the rim in the 



course of dewettin g ti me as observed by Reiter, Damman 
and coworkers [i, M [13 (for simple and for viscoelas- 
tic liquids one would expect a monotonous increase of 
the size of the rim). This feature could be understood 
qualitatively using scaling arguments and modeled by as- 
suming the friction with the substrate to be nonlinear 
in the velocity (motivated by the polymer-polymer in- 
terface at the substrate) and by adding residual stresses 
inside the film that slowly relax during the dewetting 
process [H, [13, [l3|. These stresses are argued to orig- 
inate from the spin-coating process, where upon evapo- 
ration of the solvent polymers are 'frozen in' into out-of- 
equilibrium configurations which induce an internal lat- 
eral stress [H, [13, [13] ■ When the film is heated above 
the glass temperature and starts to dewet, these residual 
stresses are able to relax, at least partially, and hence 
infiuence the dewetting dynamics. The existence of such 
stresses (for rather thick films where the occurrence of 
a shift of the glass transition temperature could be ex- 
cluded) has been shown unambiguously [l^ and they 
have been estimated to be of the order of 10^ Pa. The 
investigation of such stress-driven dewetting dynamics is 
not only of importance for the stability, the mechanical 
properties and the dynamics of thin films, but has even 
broader implications for the understanding of confined 
materials and their glass transition [U, [13]. Indeed, the 
dynamics of the dewetting has been shown to be severely 
influenced by the a ge o f the sample, both on solid and 
lic^uid substrate [l^ Il6l]. 

While there is nowadays a rather good understanding 
of the instability mechanisms leading to the film rupture 



23 , [24] , of the rim morphologies in the mature regime 
^25 , and of the thin film equations for viscoelastic fluids 

IT 11, [13, [11] , the temporal evolution of the dewetting 
process is still rather unexplored. In this work we inves- 
tigate numerically the model introduced by Vilmin et al. 
[l3j with focus on the stress-driven dewetting dynamics. 
Special interest is led on how the non-monotonous be- 
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havior of the width of the rim arises - which could not 
be explained satisfactorily so far by simple scaling laws 
[isl . [2^ which assume that friction is the only relevant 
dissipation mechanism and that the film profiles during 
dewetting are self-afhne. Therefore, we also study the 
temporal evolution of the profiles, where the viscoelastic- 
ity shows to have important effects, as well as the mutual 
importance of the two dissipation mechanisms - dissipa- 
tion by friction with the substrate and viscous dissipation 
inside the film - in the course of time. 

The model under investigation has been kept simple 
and focuses on three main features: friction on slippery 
substrate, viscoelasticity and residual stress. It corre- 
sponds to the strong slip lubrication model of Ref. [s^l in 
the limit of zero Reynolds number and without Laplace 
pressure. Since we neglect the latter, the model only 
applies to the time window before the round-up of the 
'mature rim' is appreciable, which for the highly viscous 
polymer films under consideration takes several reptation 
times. Also, instead of a hole geometry, we focus here 
on the simpler case of an edge geometry, i.e. a straight 
contact line. Then the system can be described by a 
one-dimensional model. Experimentally, the dynamics 
in the hole and the edge geometry have been compared 
in Refs. 0, [3l| and a generalization of the model to the 
hole geometry is possible, see Ref. [29^. 

This work is organized as follows: In section [TTl we 
briefly review the model and the assumptions made in 
its derivation. In section lllli we give analytical expres- 
sions for the short-time behavior of the model, both as 
a benchmark for the subsequent numerical work and to 
introduce the characteristic length and velocity scales. 
Section IIVI reformulates the problem in dimensionless 
variables, the reduced parameters are discussed and the 
numerical method is briefly described. Section |V] shows 
numerical results for the physical observables: the dewet- 
ted length, the velocity at the edge, the height of the rim 
at the edge and the width of the rim. In the following 
section IVll we have a closer look at the height, velocity 
and stress profiles. Section fVIII is devoted to the energy 
balance in the course of dewetting that is used to extract 
the two relevant dissipation mechanisms from the nu- 
merical results. Finally we discuss the interplay between 
viscoelasticity, friction and residual stress in section IVIIII 
and conclusions are presented in section [IXl 



II. MODEL 

We consider dewetting from a straight edge to get an 
effective one-dimensional description. Fig. [T] shows a 
sketch of this geometry. Initially the film is flat with 
a vertical front at x = and the film extends infinitely 
in the direction x > 0. Since we are interested only in 
the dynamics of the rim for times smaller than the time 
scale of the build-up of the mature rim, we neglect the 
Laplace pressure arising from the film-air interface cur- 
vature. The film thickness t) is assumed to be small 




x=0 



x=L(t) 



FIG. 1: Sketch of the thin film geometry. The height profile, 
h{x, t), is represented in grey. The dewetted length, L{t), and 
the width of the rim, W{t), are defined. The velocity, v{x,t), 
is shown as black arrows; it decreases with the distance from 
the edge. The residual stress, o-o(t), is indicated as a white 
arrow; it is assumed homogeneous in space, but relaxes in 
time (see text). 



with respect to the hydrodynamic extrapolation length, 
or slippage length, b = 77/ where rj is the viscosity and 
C is the friction coefficient with the substrate [l^. This 
is the situation of interest for polymer films on polymer- 
covered substrates and results in a plug flow in the film. 
In the spirit of a lubrication approximation, it is then 
sufficient to consider the horizontal velocity v{x,t), i.e. 
in the direction perpendicular to the dewetting front. For 
details concerning the derivation and the lubrication ap- 
proximation we refer to Ref. [29j . 

In brief, the horizontal momentum equation (inte- 
grated over the thickness of the film) is the balance of 
the frictional force at the film-substrate interface and the 
divergence of the total stress inside the film. We allow 
for cither linear friction, — (^v, or nonlinear friction, 



(1) 



as motivated by recent experiments on polymer-polymer 
friction [s^l- Here ti is a characteristic velocity above 
which the friction is nonlinear and a is the exponent 
characterizing the nonlinear behavior of the friction. If 
not specified otherwise we use a = 0.8, a value ob- 
tained recently from experimental data on PS dewetting 
on PDMS-covered substrates [l3|, which is also in qual- 
itative agreement with measurements on rubber-brush 
friction [33 • The divergence of the stresses in the film 



is F, 



_ d(ha) 
dx 



Inertia is neglected due to the high vis- 
cosity of the polymer film, implying overdamped motion. 
Momentum balance then leads to 



djha) 
dx 



(2) 



with a = for linear friction and, in general, < a < 1 
for nonlinear friction. 

To connect stresses and velocity gradients, a consti- 
tutive relation has to be specified. We account for the 
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viscoelasticity of the polymer film by using a standard 
Jeffrey model I33l |. 



■ + Tidt<T = Gti [{d^v) + Todt (d^v)] 



(3) 



Here G is the elastic modulus and tq and ti , with tq <C ti , 
are the two characteristic time scales. By use of these 
one can define two viscosities r?o = tqG = 2t] and 
771 — TiG. The constitutive relation Eq. ([3]) describes vis- 
cous behavior with viscosity 770 for times t < tq, elastic 
behavior for times tq < t < ti and again viscous behavior 
with viscosity 771 for times t > ti. The latter, highly vis- 
cous behavior (since rji ^ rjo) at long times is supposed 
to reflect the slow dynamics by polymer reptation. The 
time scale ti may be substantially smaller than the bulk 
relaxation time, due to the thin film geometry .35| . but 
has been recently proven to scale with molecular weight 
like the reptation time |17i] . The short-time dynamics, 
t < tqi involves intra-chain motions and is dominated by 
monomer friction resulting in a much lower viscosity. 

The change in the height profile of the film is given by 
volume conservation 



dth = -dx (vh) . 



(4) 



The position of the dewetting front is L{t), starting from 
L{t = 0) = 0. Its motion is governed by the velocity in 
the film at the edge, i.e. by dtL — v{L). 

Now we have to specify boundary conditions. At the 
edge of the film, the height-integrated stress has to equal 
the driving force, hence 



h{L)a{L)^~\S\. 



(5) 



Here S = jsv is the so-called spreading parameter 
and 7s«, 7sz and 7 are the surface energies for solid- vapor, 
solid-liquid and liquid-vapor interfaces respectively. In 
our case we consider a non-wettable substrate, i.e. < 0. 
Additionally we impose v{x = 00) — 0, assuming that 
the dewetting front is far away from other fronts or holes 
and that the film is unperturbed far away. As initial 
conditions we prescribe v{x,t = 0) = and h{x,t = 
0) = /iQ, i.e. a quiescent film of height /iq. The initial 
stress a{x,t = 0) is either zero, which we refer to as the 
case without residual stress, or given by a constant value 
throughout the film, cto, which we refer to as the case 
with residual stress. 



III. ANALYTICAL SHORT-TIME SOLUTION 

One now has to solve Eqs. (OH]) with the initial and 
boundary conditions specified above. In general, this has 
to be done numerically. For short times however, one 
can give insightful analytical expressions, due to the facts 
that the height profile is initially only slightly perturbed, 
\h{x) — ho\ <C ho, and that the system is purely viscous 
at very short times, i.e. a ~ 2rjdxV. 



A. Without residual stress 

We first focus on the case without initial residual 
stress, i.e. with a{x, i = 0) = 0. To leading order one 
has to solve 



(6) 



which in the case of linear friction, a = 0, yields an 
exponential velocity profile 

v(x, t) — Vn exp I — ^ I for x>L. (7) 

V V2W0J 

Here L = L{t) denotes the dewetted distance, with Lit = 
0) = 0. Wo is the characteristic rim width. By matching 
the boundary condition, Eq. (O, one can identify the 
characteristic initial velocity Vq . These two characteristic 
scales read 



Vo 



Wo 




(8) 



The characteristic rim width can be rewritten as Wo = 
^/hob, where b is the extrapolation or slippage length, an 
expression that has been proposed some time ago |36l.l37l|. 

In the case of nonlinear friction, a ^ 0, a, velocity 
profile of the following form is obtained: 



/ a X — L 

V[X,t) = Vo,a 1-77- 



2 V2Wo,a 



(9) 



for < a; - L < ^V^Wo^a and v{x) = elsewhere. The 
corresponding characteristic scales now read 



2 \2 



Wo 



V' 



2-a\ Vq" 



and 



(10) 



with Vo and Wq as defined in Eq. ([8]) above. 

For the height profile at short times one simply has to 
solve dth = —hodxv{x) which results in 



h{x,t) — ho 
h{x,t) = ho 



1 



exp 



^/2Wo 

\/2Wo,a 



X — L 
'V2W0 



and 



1 - 



a X 



2 V2W0.C. 



(11) 



for linear and nonlinear friction, respectively. In both 
cases, the rim build-up is initially linear with time, 
h{L{t),t) oc t. The stress inside the film is given by 
(j{x) — 2rjdxv{x) and it can be easily verified that at the 
edge hocr{L) = —\S\ holds, and that the stress vanishes 
for X 00. 
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Effect of residual stress 



Thus we arrive at the dimensionless equations 



We now consider the case of an initial residual stress, 
a{x,t = 0) = (To- Since the constitutive law is linear the 
total stress can be written as a*^°^{x) — <t^{x) + CTq (see 
sections |V] and IVIIBI for a numerical validation) . Here 
cr^ represents the stresses in the absence of the residual 
stress (the superscript v stands for viscoelastic) . For the 
short-time response, we can again use = 2r]dxV and 
the stress balance at the edge reads 



'\L 



= 2ridxv{x)\i 



(To = 



1^1 
h{L) 



(12) 



Far away from the edge, there are no flows and addition- 
ally a{x = oo) = ctq must hold. Using these boundary 
conditions, one gets the same formulas for the velocity 
and height profiles, i.e. Eqs. jT]), ^ and Eqs. HI]) for 
linear and nonlinear friction respectively, but with the 
substitution 



\S\^\S\+h„a„ 



(13) 



entering the characteristic velocity scale Vq, and in the 
case of nonlinear friction entering both Vq q, and Wq^q. 

This clearly indicates that, as expected, the resid- 
ual stress constitutes an additional driving force for the 
dewetting. For both linear and nonlinear friction, in the 
presence of a residual stress ctq the initial velocity of the 
dewetting process is increased by a factor {1 + hoao/\S\) 
and (1 -I- ho<To/\S\)^^ , respectively. The initial stress 
profiles in the presence of residual stress read 



a{x) 
a{x) 



(^0 




for a = and < a < 1, respectively. 



IV. DIMENSIONLESS EQUATIONS, 
PARAMETERS AND NUMERICAL METHOD 

To get access to the dynamics at longer times, one has 
to solve Eqs. numerically. For this purpose, it is 

convenient to rescale these equations. Keeping in mind 
the typical length and velocity scales obtained above, we 
define the following dimensionless quantities 



h'^^,t'^-. (15) 
ho To 



Wo,c. ' 



At short times, a = 2r]dxV holds, so we rescale the stress 



as a' = a / a* with a* = 2rj 



Vq, 
Wo,, 



2dx' {h'a') , 



(16) 



a' + T[dt'<j' 
dt'h' 



T[[{dx'v')+dt'{dx'v')] , (17) 

-pdx' iv'h') , (18) 



with t[ = ti/tq. 
parameter 



We also have introduced an reduced 



Wq 



= V2 



Gho 



(19) 



which is a dimensionless number quantifying the coupling 
strength of the flow field to the height field, and which is 
proportional to |S'|/G. 

The rescaled boundary and initial conditions read 



1 



h\L)a'{L) = - — 



(20) 



v'{x = oo) = 0, v'{x',t' = 0) = 0, h'{x',t' 0) = 1 
and a'{x',t' = 0) = or ctq respectively. For the sake of 
simplicity, in what follows we will suppress the primes. 

Remarkably, we are left with only three effective pa- 
rameters: first, the friction exponent a that governs 
whether one deals with linear or nonlinear friction. Sec- 
ond, the time scale r( (in units of the short time scale 
To) which governs the viscoelastic crossover from the elas- 
tic behavior of the film towards the highly viscous flow 
regime by reptation of chains in the film. Finally, the 
parameter f3, describing the coupling between flow and 
height profiles. 

We solved Eqs. (fTB)) - P^ numerically, with the bound- 
ary and initial conditions specified above. For this pur- 
pose, we discretized time and space. At every time step, 
we applied a shooting method to solve Eqs. ([16]) and (fTT)) 
simultaneously: i.e. at the edge we started with a trial 
velocity v*{L) (e.g. the one of the last time step) and 
the stress value prescribed by the boundary condition, 
Eq. (PO]) . and evolved both equations in space. This was 
done iteratively adjusting v*{L), until the boundary con- 
dition v{x^l) = was satisfied up to a prescribed tol- 
erance. Then we updated the height profile on a moving 
grid and proceeded to the next time step. 

In this work we are not so much interested in the ef- 
fects of parameter variations, but instead want to focus 
on generic properties and the interplay between viscoelas- 
ticity, linear or nonlinear friction, and the relaxation of 
residual stress. Hence we use throughout the paper the 
parameter values ri = 100, (3 = 0.2 and a = or a = 0.8 
for linear and nonlinear friction, respectively. The initial 
residual stress was set either to a{t = 0) = or 1. 

A difference of two orders of magnitude in the two time 
scales To and ti is enough to separate the two viscous 
regimes and is a good compromise to avoid too time- 
consuming simulations. The actual value of ti/to is yet 
unknown because the short time scale To is hard to ac- 
cess experimentally. The value for the nonlinear friction 
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FIG. 2: (Color online) The dewetted length, L{t), is shown 
in double logarithmic scale. Nonlinear friction with a = 0.8 
was used. The black solid line shows the dependence in the 
absence of initial residual stress, while the red dashed curve 
has been obtained for a{t = 0) = ao = 1. The inset shows 
the curves in linear scale. 

exponent, a — 0.8, was used following the experimen- 
tal results from Ref . fisj . The value of the parameter (3 
has been chosen for numerical convenience. This param- 
eter has only quantitative effects on the dynamics, and 
modest variations are noticeable only in the evolution of 
the height profile. However, we use reduced variables - 
the elastic modulus G and the wetting parameter \ S\ that 
enter /? of course do have an important influence: the for- 
mer is connected to both viscosities, and thus both enter 
the scales of the velocity and the width, see Eqs. (|8ll0p . 



V. NUMERICAL RESULTS FOR PHYSICAL 
OBSERVABLES 

In dewetting experiments, two observables are directly 
accessible by optical microscopy |9i] : the dewetted length 
L{t) and the width of the rim W{t). The velocity at the 
edge is then obtained by taking the temporal derivative 
of the dewetted length, ^ = Vit) = v{L{t),t). The 
height of the film at the rim is not that easily accessible 
(usually atomic force microscopy has to be used) and a 
full picture of the temporal evolution is hard to obtain. 

In our simulations, we have direct access to these ob- 
servables. Fig. [2] shows the dewetted length L{t) in a 
double logarithmic scale (see the inset for linear scale). 
Fig. [3] displays the temporal evolution of the velocity at 
the edge, V{t) = w(L), again in double logarithmic scale. 
In both figures, the black solid lines are the results with- 
out initial residual stress. Since we have rescaled the ve- 
locity by Vb,a, y{i) starts at one. Once the elastic regime 
is entered (for t > tq = 1) the velocity decreases rapidly: 
the behavior is roughly t^^ (see the dotted line in Fig. [3]), 
as is discussed in detail in Ref. [l^. Upon reaching the 
long time viscous regime (t > ti = 100), the velocity still 
decreases, but only slowly. The same behavior can be de- 
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FIG. 3: (Color online) The velocity at the edge, Vit) = w(L), 
is shown as a function of time in double logarithmic scale. 
Corresponding to Fig. [2j the black solid line is in the absence 
of residual stress and the red dashed curve has been obtained 
for cr(t = 0) = ao = 1. The green dash-dotted curve shows 
a control simulation with a residual stress that is not allowed 
to relax (see text for details). 

duced from Fig. [21 where the dewetted length is almost 
linear in time, actually slowly decreasing, for t < 1 and 
t > 100. 

The red dashed lines in both figures have been obtained 
with an initial residual stress, a{t = 0) = (Tq = 1- As 
expected, since the stress acts as an additional driving 
force (see Eq. ([T51) ). the initial velocity is much higher 
in the presence of the residual stress. However, it also 
decays more rapidly due to the relaxation of this extra 
driving, until after a time of order ti both velocities are 
approximately the same. 

Figure 2] shows the temporal evolution of the height 
of the film at the edge, H{t) = h{L{t),t). The black 
(solid) and red (dashed) curves are obtained without ini- 
tial stress and with a{t = 0) = 1 and correspond to 
the respective velocities displayed in Fig. [31 One can 
clearly discern the two viscous regimes at short and long 
times and the intermediate 'elastic plateau'. The resid- 
ual stress leads to more pronounced rims, i.e. higher 
values of h{L). This could be expected from the short- 
time analysis, Eqs. pT|) and ifTS)) . It remains true for all 
times, because there is no mechanism that could lead to 
a decrease of the height, even if the residual stress has 
relaxed. This figure remains unchanged if linear friction 
(a = 0) is used instead. 

Figure [51 shows the temporal evolution of the width of 
the rim, W{t). It has been obtained by introducing a 
cutoff as proposed previously [l^: W{t) is defined as the 
distance where h{x,t) = (1 -I- l/10)/io. For usual liquids 
one would expect a monotonous increase in the width, 
even in the mature regime. This is the case without ini- 
tial residual stress (black solid curve) as well as with an 
initial stress that does not relax (green dash-dotted curve, 
see below). Only if the residual stress relaxes, here with 
the time scale ri given by the Jeffrey model, in the course 
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FIG. 4: (Color online) The height at the edge, H{t) = h{L), 
is shown as a function of time in double logarithmic scale. 
The black (solid) and red (dashed) curves correspond to the 
respective curves displayed in the previous figures. They are 
obtained with a{t = 0) = and 1, respectively. One can 
clearly perceive the elastic plateau. 



of the dewetting process a maximum in the rim width is 
obtained as shown by the red dashed curve. This max- 
imum is located close to the time scale ri. The value 
of the maximum increases with increasing residual stress 
tJoj as has been studied previously ^2§\. 

We note that in this work we implement the residual 
stress as an initial bulk stress, a{t = 0) = ctq, see sec- 
tion ini This implies that the dynamics of this stress is 
governed - as are the viscoelastic stresses - by the consti- 
tutive law, Eq. pT|) . and hence it decays with the time 
scale Ti. The green dash-dotted lines in Figs. [3] and [H 
have been obtained by a simulation where the residual 
stress was implemented differently: As has been already 
proposed in section lTlI Bl since the constitutive law is lin- 
ear one can separate the two stresses right from the start 
by writing a^°*-{x) — cr^ix) + ctq. The residual stress, 
i.e. the contribution ctq, then appears as the additional 
driving as described by Eq. and additionally as a 

homogeneous contribution in the momentum equation, 
Eq. (|16p . One can use this separation, impose zero initial 
bulk stress and consider ctq as a parameter that explicitly 
decays like ao{t) = ao{t = 0)e~*/'^i. Indeed, in doing so 
one exactly regains the results obtained by implement- 
ing an initial value for the stress in the bulk, i.e. the 
red dashed curves, confirming that this separation in the 
stress is consistent. In contrast, the green dash-dotted 
curves in Figs. [3] and [S] have been obtained by imposing 
fJo = 1 for all t > 0. This represents the case where 
the residual stress does not relax (or only on time scales 
much larger than the time scale ti). Initially, the veloc- 
ity at the edge and the width of the rim are the same 
for this case and for relaxing stresses, see the green and 
red curves, hence again confirming the additional driving 
by the stress. However, there is no maximum in the rim 
width but a monotonous increase if the residual stress 
does not relax. 



W(t) - 
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FIG. 5: (Color online) The width of the rim, W{t), as a func- 
tion of log(t) for the velocity evolutions shown in Fig. [3] The 
black solid line, obtained without initial residual stress, and 
the green dash-dotted line, obtained with a non-relaxing resid- 
ual stress (see text for details), display a monotonous growth 
of the width of the rim. Only the red dashed curve, obtained 
with a residual stress that relaxes (here with time scale ri) 
has a maximum that is located close to n . 

At this point we can conclude that the model with 
residual stress reproduces most of the features of the 
experiments in Refs. [1, [l^ , namely the high initial 
velocity, the fast decay of the dewetting velocity approx- 
imately like v{L) ~ t~^, the elastic plateau in the height 
of the rim and the maximum in the rim width. It has 
been shown to be crucial that this residual stress relaxes. 
This relaxation dynamics then naturally influences the 
dynamics of the film. To get a better understanding of 
the processes involved especially in the nonmonotonous 
behavior in the rim width, in the following sections we 
focus on the temporal evolution of the profiles and the 
dissipation mechanisms involved. 



VI. NUMERICAL RESULTS FOR THE 
PROFILES 

For the system we are modeling, PS on PDMS-covered 
substrate, the evolution of the height profiles of dewetting 
films has been studied in the short-time regime by atomic 
force microscopy Q, revealing very asymmetric profiles. 
In principle, the height profile could be extracted also by 
investigating the interference patterns in optical micro- 
graphs. Although velocity fields have been made visible 
e.g. in bursting of suspended soap films [38], the velocity 
profiles are probably not obtainable by simple means for 
a highly viscoelastic polymer film. 

In our simulations, we have direct access to the profiles 
and their temporal evolution. Fig. [S] shows the height 
profiles of the film for successive times. The upper panel 
displays the case without residual stress, while the lower 
panel was obtained with a{t = 0) = 1. At short times, 
we get indeed the profiles as described by Eq. (fTT|) , in the 
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h(x) 




FIG. 6: (Color online) The height profiles h{x) at successive 
times: black (solid) curve t = 0.1, red (dashed) t = 1, green 
(dotted) t = 10, blue (dash-dotted) t = 50, orange (long 
dashes) t — 200. The upper panel shows the case without 
residual stress, the lower panel has been obtained with an 
initial residual stress a{t = 0) — 1. 



case with residual stress with the substitution Eq. (fT3|) . 
At longer times, the form of the profiles can not be given 
by a simple function anymore. In accordance with Fig.[5l 
in the case without residual stress the width of the profile 
is steadily increasing, while in the presence of the residual 
stress, there is a regime where the edge advances faster 
than the other end of the rim and the width changes 
non-monotonously. 

In the course of the dewetting process, the velocity pro- 
files change more severely than the height profiles. Fig.[7| 
shows the velocity profiles inside the film for successive 
times in the absence of initial residual stress, and Fig. [5] 
for the case of a{t = 0) = 1. In both figures the veloci- 
ties have been renormalized to 1 at the edge, in order to 
make the changes in the shape visible. The insets show 
the unrenormalized velocities that are rapidly decreasing 
in amplitude, as can be seen from Fig. [3] for the (max- 
imum) value at the edge. The short-time profiles are 
in accordance with the predictions of Eqs. (O and 
The velocity as a function of the distance from the edge, 
X — L, is rapidly decaying. However, when the system 
is in the elastic regime (tq < t < ri), the profiles are 
severely perturbed and become concave, see the green 
dotted curves in both Fig. [7] and Fig.[51 After the system 
has left the elastic regime and has again become viscous, 
for t > Ti, the behavior is distinct for the two cases, see 
the blue (dash-dotted) and orange (long-dashed) curves: 
without residual stress the shape of the profile becomes 
convex again with similar (actually still slowly increas- 
ing) widths. With residual stress, the profiles also regain 
a convex shape, but first with a large characteristic width 
that then retracts. 

The concave, rounded shape of the velocity profiles in 
the elastic regime can be understood qualitatively from 
the governing equations. In the elastic regime, Eq. ([3]) 



FIG. 7: (Color online) The velocity profiles v{x), without 
residual stress, at the successive times as in Fig. |S] The ve- 
locities at the edge have been normalized to one to make the 
changes in the profile visible. The inset shows the actual pro- 
files, where the edge velocity is rapidly decaying, cf. Fig. [S] 



v(x) 




FIG. 8: (Color online) The velocity profiles v{x), with an 
initial residual stress of a(t = 0) = 1. The successive times 
are again as in the two previous figures and the inset shows 
the unrenormalized, rapidly decaying velocities. 

states that the dominating contribution to the stress is 
dtu — GdxV. If we assume linear friction for simplicity, 
for the velocity we have Ci' = dx (ha), leading to 

dtv^C'Gdxihpdxv). (21) 

Here we have used the fact that the height is almost time- 
independent in the elastic plateau, see Fig. [H and has an 
almost stationary profile (in the frame moving with L) 
that we call hp{x). Eq. (PT|) is a generalized diffusion 
equation for the velocity field. The initial condition to 
think of is (upon entering the elastic regime) an exponen- 
tial profile, cf. Eq. ([7]). Additionally, there is a boundary 
condition at the edge, namely that v{L) is decaying in 
time, cf. Fig. [31 It is easy to convince oneself that in- 
deed the exponential evolves towards a concave, rounded 
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o(x) 




FIG. 9: (Color online) The stress profiles cr{x), with an initial 
residual stress of a{t = 0) = 1. The successive times are: 
black (solid) t = 0.1, red (dashed) t = 1, green (dotted) t — 
10, blue (dash-dotted) t = 50, orange (long dashes) t = 500, 
violet (dash-double dot) t = 1000. 



right hand side has to be so. Integration by parts of the 
right hand side, making use of the boundary conditions 
which imply [hav]'^ — —h{L)a{L)v{L) = \/2v{L), and 
rearrangement of terms results in 

V2v{L)^— v^-"dx + 2 had^vdx. (23) 
2 Jl Jl 

This formula has a simple interpretation: the left hand 
side is the work done by the rescaled driving force. In 
dimensional form it is proportional to |iS'|w(L) and thus 
has units of force/length times velocity. It naturally ap- 
peared as a boundary term. On the right hand side, the 
first term is the dissipation by the (in general nonlin- 
ear) friction of the film with the substrate, and the sec- 
ond term is the height averaged (hence the weight factor 
h{x)) dissipation inside the film, which is of the usual 
form adxV. 



A. Effect of residual stress 



shape in the course of time. This happens the more 
severely the larger G/Q is, and the faster v{L) decays 
- which is the reason for the shape being more severely 
perturbed in the case of residual stress, where v{L) de- 
cays much faster, see Figs. [7] and HI The spatially de- 
pendent 'diffusion coefficient' hp(x), that decreases with 
the distance from the edge, is enhancing the rounding up 
close to the edge. The same general argument holds in 
case of nonlinear friction, although it is less obvious since 
Eq. (I2ip becomes a nonlinear diffusion equation. 

The profiles for the stress inside the film are shown in 
Fig. [51 in the presence of residual stress. The short time 
profile is in accordance with Eq. (|14p . At the edge, the 
stress is determined by the boundary condition, Eq. (j20p . 
The relaxation of the initial residual stress can be de- 
duced from the value at large x. Relaxation is complete 
after several relaxation times ri . 



VII. ENERGY BALANCE 

Since the work described in Rcfs. [ll,!!^!, it has been 
proven successful to describe the temporal behavior of 
dewetting by an energy balance that accounts for the rel- 
evant driving and dissipation mechanisms. This method 
has also been recently used to study the dewetting in the 
present model by means of scaling laws [29j . The energy 
balance can be derived from Eq. (jl6p by multiplication by 
v{x) and subsequent integration over space, which yields 



In the presence of residual stress, we again write 
a*-°^{x) = a-''{x) + (To. From Eq. we thus obtain 

v'--dx 



22) 



v{xf~°'dx^2 / dx(^h{x)a{x)^v{x)dx. ( 



For brevity we have suppressed the temporal dependence 
in the fields v, h and a, but this equation has to hold for 
all times t > 0. The left hand side is the rescaled dissi- 
pation by friction. It is obviously positive, thus also the 



-V2v{L) + 2 ha^'dxvdx + 2ao hd^vdx . (24) 
Jl Jl 

where we used that the residual stress is considered ho- 
mogeneous here. Both the second term (proportional to 
~v{L) < 0) and the last term are negative (since in gen- 
eral dxV < 0, while h > and ctq > 0). Thus these 
terms should be interpreted as the work done by the in- 
creased driving force due to the residual stress [33]. The 
remaining terms are positive (for the third term one has 
to note that CTv < and dxV < 0) and as before can be 
interpreted as the dissipation by friction and the viscous 
dissipation inside the film, respectively. In total, we get 
the balance 

V2v{L) — 2cto / hdxvdx 
~" ' v'^-°'dx + 2 I ha^dxvdx. (25) 



It is of the same form as Eq. ([23]) , but with the additional 
driving arising from the residual stress. A similar relation 
was introduced empirically in Ref. [HI. Note that since 
the dynamics of the residual stress is governed by the 
constitutive relation, Eq. ITTI) . ctq in Eq. ([25)) is time- 
dependent and relaxes like cro(i) = (To(O)e^*' '^^ . 



B. Numerical evaluation 

It has been shown for viscous fluids [s^] that if a rim 
has already been formed on a slippery substrate, the vis- 
cous dissipation should be negligible as compared to the 
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FIG. 10: (Color online) Numerical evaluations of the energy 
balance, Eq. (|25|l . The black solid curves show the work done 
by the driving force, AVI^, the red dash-dotted curves the 
dissipation in the film, TSv, and the blue dashed curves the 
dissipation by friction with the substrate, TSf. Panels (a) and 

(b) correspond to the case with linear and nonlinear friction 
law, respectively, in the absence of residual stress, while panel 

(c) shows the case of nonlinear friction and an initial residual 
stress a{t = 0) = 1. 



dissipation by friction with the substrate. Recently it has 
been shown |13j . that for a linear friction law and in the 
initial (viscous) regime of dewetting, both mechanisms 
contribute equally to the dissipation. What happens at 
intermediate times, in the viscoelastic regime, and what 
are the effects of nonlinear friction and residual stress is 
explored in this section. 

Knowing all the relevant fields from the numerical solu- 
tion, using Eq. (or Eq. ([^5]) in the presence of resid- 
ual stress) we have access to the work done by the driving 
force, AW = ^/2v{L) — 2aoJpidxvdx, the dissipation by 
friction TSf — J^^~"dx and the dissipation inside 

the film TSy — 2 J^a^dxvdx. These three quantities 
are shown in Fig. [TO] the upper panel (a) shows the case 
with linear friction {a = 0) and without residual stress. 
The middle (b) and lower (c) panel have been obtained 
with nonlinear friction, without residual stress and with 
a{t = 0) = 1, respectively. In all panels the black (solid), 
red (dash-dotted) and blue (dashed) lines correspond to 
AW, TSv and TSf respectively. As it should, the sum 
of the blue and the red curves, i.e. the total dissipa- 
tion, equals the black curve, which is the work done by 



the driving force, confirming that the energy balance is 
respected. 

At short times and in the case of a linear friction law, 
Fig. [TOT a) shows that the two dissipation mechanisms 
contribute equally, in accordance with Ref. [l3l |. This 
can also be seen by directly evaluating Eq. (^5)) with the 
analytical short-time solutions obtained in section HTTl In 
the elastic regime, for tq < t < ti, the dissipation by fric- 
tion is larger than the dissipation in the film. Since the 
system is rather elastic than viscous, the dissipation in 
the film is reduced, until it becomes again comparable 
to the frictional dissipation for t > ti. Thus the dis- 
sipation by friction is always more or at least equally 
important than the dissipation inside the film. This is in 
contrast to the case with nonlinear friction displayed in 
Fig. fTOT b) and (c). In this case, at short times the viscous 
dissipation is dominating, since the sublinear friction law 
reduces the friction. In the elastic regime, the frictional 
dissipation becomes dominating, by the same argument 
than in case of linear friction. For t > ti, where the 
system becomes viscous again, the dissipation in the film 
is dominating, especially in the case with residual stress 
shown in Fig.fTUTc). This change in the importance of the 
two dissipation mechanisms is clearly a combined effect 
of both viscoelasticity and nonlinear friction. 



VIII. DISCUSSION 

We are now able to discuss the effects of the viscoelas- 
ticity, the nonlinear friction and the initial residual stress 
and their mutual interplay in some detail. The viscoelas- 
ticity clearly is responsible for the plateau in the height of 
the rim, cf. Fig. [4l that is also experimentally observed. 
Concerning the flow profiles, it leads to a rounded con- 
cave form of the velocity profile, cf. Fig. [7] that can be 
explained by the elasticity and the decrease in the driving 
force due to the rim build-up. When the system enters 
the second viscous regime, around ri, the elastic energy 
stored in the fiow profiles is dissipated in the film, leading 
to a relative increase of the viscous dissipation with re- 
spect to frictional dissipation. In case of a linear friction 
law with the substrate both dissipation mechanisms, the 
one by friction on the substrate which was dominating 
in the elastic regime and the one by viscous dissipation 
inside the film, become approximately equally important 
for t > Ti, cf. Fig.fTUTa). In case of nonlinear friction the 
dissipation in the film becomes even more important than 
the dissipation by friction, see Fig. [Torb').(c). Though, 
on even longer time scales (not captured by the model 
presented here) where the mature regime is reached and 
the profile becomes a half-cylinder with a stationary flow 
profile while the rim is still growing, the dissipation by 
friction will be dominating again ^37]. 

The nonlinear friction renders all the profiles steeper, 
i.e. those for the height of the film as well as those 
for the velocity and the stress inside the film. Together 
with the viscoelasticity, the nonlinear friction leads to 
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two crossovers concerning the importance of the two dis- 
sipation mechanisms: while in the two viscous regimes 
(for t < tq and i > n) the viscous dissipation in the fihn 
is more important, since the friction is reduced due to 
the subhnear friction law, in between tq and ri the dis- 
sipation by friction on the substrate is dominating. This 
is due to the fact that the polymer film is predominantly 
elastic in this regime and viscous dissipation is reduced. 

The residual stress leads, as expected, to an increase of 
the driving force, cf . Eq. ([T3|) , consequently resulting in a 
faster initial dynamics. Concerning the velocity profiles, 
the faster decrease of the driving force due to the relax- 
ation of the initial stress amplifies the roundup of the pro- 
files, cf. Figs. [7] and m Hence, when the system changes 
from elastic to viscous around ti, the dissipation inside 
the film is larger than without stress. This indicates that 
the occurrence of the maximum in the rim width is due 
to the faster decrease in the driving by stress relaxation, 
coupled to the increased importance of dissipation inside 
the film that was caused by both the viscoelasticity and 
the nonlinearity in the friction. 

IX. CONCLUSIONS AND PERSPECTIVE 

In conclusion, we have investigated numerically the 
temporal evolution of a dewetting thin polymer film un- 
der homogenous lateral stress. The main properties of 
the polymeric system, namely the viscoelasticity of the 
film, the friction with the substrate and the residual 
stress have been shown to interplay in the course of 
dewetted time, leading to complex behavior. Numeri- 
cal evaluations of the energy balance, derived from the 
underlying model equations, revealed that the dissipa- 
tion during the dewetting is more complex than expected. 
The viscous dissipation in the film and the frictional dis- 
sipation at the substrate have different time dependences 
and, in case of a nonlinear friction law, the friction is not 
always the most important dissipation mechanism. Es- 
pecially, the occurrence of a maximum in the rim width, 
as observed experimentally in Refs. [IB, [13]: could be 
traced back to the combined effect of a more rapid de- 
crease in the driving force due to relaxation of stress 
and the viscous dissipation of the elastic energy stored in 
the flow profile, which is the dissipation mechanism with 
the largest contribution at the time scale of the elastic- 
viscous crossover due to the nonlinearity of the friction. 

Our numerical treatment also sheds new light on sim- 
ple approaches based on scaling laws put forward pre- 



viously [H, [2^, [131 ■ There one either assumes that the 
dissipation by friction dominates over the viscous dissipa- 
tion, or that both mechanisms have the same temporal 
dependence. Additionally, in these approaches one has 
to assume a simple form of volume conservation, namely 
hoL{t) = C{H{t) — hi))W{t), to connect the rim width 
W with the dewetted length L and the height H at the 
edge, with a fixed constant C depending on the shape of 
the rim. The complex behavior concerning the dissipa- 
tion and the non self-similarity of the profiles indicates 
that these simple scaling arguments should be revisited. 
Indeed, if dissipation by friction would be the only im- 
portant dissipation mechanism, by balancing the driving 
with the frictional dissipation, a decrease in the width 
of the rim would result in a speed up of the dewetting 
velocity, which is never observed. 

The inclusion of a residual stress in the model proved 
to be crucial for the occurrence of the maximum in the 
rim width. This maximum has been already used to ex- 
tract relevant information from experiments, namely the 
exponent of the nonlinear friction law [l8| . The inclusion 
of residual stresses has also been put forward to analyze 
effects of film ageing: The dynamics of the dewetting has 
been shown to be severely infiuenced by the age of the 
sample, both on solid and liquid substrate [l3. Il6j. 

Clearly it would be desirable to extract more informa- 
tion on these stresses both from the technological point 
of view of film stability and for fundamental reasons to 
better understand the ageing of a confined glassy poly- 
mer film by means of dewetting experiments. Less is 
known so far on how these stresses look like - what are 
the non-equilibrium configurations and the relaxation dy- 
namics of the polymer chains confined in such a thin 
film? There is no special reason, a priori, for the resid- 
ual stress to relax with the same time scale that governs 
the long-time flow behavior of the polymeric liquid - as 
we have assumed here for simplicity. Indeed, it has been 
recently observed experimentally ^17] that there are (at 
least) two time scales: the characteristic time of the max- 
imum in the rim width, which does not scale like repta- 
tion with molecular weight, and the long-time crossover 
in the dewetting velocity, which does. These issues will 
be investigated in a forthcoming publication. 
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